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AHS I K ACT 

Numerical  experiments  for  the  unsteady  transonic  fli>»  past  a symmetric  ntrtoil 
with  an  oscillating  flap  in  free  air  are  carried  out  for  a range  of 
supercritical  Mach  numbers  and  reduced  frequencies  using  two  newly  devised 
computat ional  schemes.  The  inviscid  results,  unhampered  by  the  complicating 
effects  of  wind  tunnel  wall  interference  and  shockwave  and  boundary  layer 
Interaction,  evaluate  the  extent  to  which  the  unsteady  loading  responds 
linearly  to  changes  in  flap  deflection  and  also  the  dependence  of  shock  excursion 
amplitude  and  net  unsteady  lift  and  moment  coefficients  on  oscillation 
frequency.  The  main  conclusions  are  discussed  in  light  of  recent  experimental 
1 ladings . 


II  1 N l'KODUlTlON 


A characteristic  of  transonic  unsteady  flow  with  shockwaves  is  the 
potentially  large  phase  lag  between  boundary  motion  and  induced  surface  pressure. 
Moreover,  net  force  coefficients  generally  exceed  those  in  subsonic  and 
supersonic  speed  regimes.  These  effects  tend  to  increase  the  likelihood  of 
aeroelastic  instability,  making  transonic  speeds  most  critical  for  aircraft 
flutter.  I.andahl's  classic  monograph  [ 1 ] addresses  the  basic  physical  and 
mathematical  issues,  emphasizing  tractable  linear  models. 

The  inherent  nonlinearity  of  the  problem  obviously  forbids  general 
analysis.  However,  two  somewhat  limited  numerical  capabilities  currently 
exist.  The  first  analyzes  the  unsteady  loading  Induced  by  small  amplitude 
oscillations  using  a relaxation  solution  of  the  mixed  type  time  Fourier 
transformed  disturbance  equation  which  results  from  linearization  about  a 
prescribed  steady  nonlinear  mean  flow  and  applies  to  motions  of  arbitrary 
frequency  [“]  • The  second,  somewhat  more  general,  integrates  the  governing 
small  disturbance  equation  in  time  using  an  efficient  alternating-direction- 
implicit  factorization,  but  with  a restriction  to  low-frequency  motions  j • 
Shortcomings  are  present  in  each  of  these  methods.  In  the  tirst  case  the 
nonharmonic  part  of  the  total  disturbance  flow  is  artificially  fixed  and  taken 
as  the  steady  sol  lit  ion  that  obtains  when  no  forced  oscillations  are  present: 
the  required  nonlinear  feedback  and  energy  transfer  from  harmonic  components 
to  mean  flow,  accounted  for  in  direct  time  integration  methods,  which  depend 
on  both  oscillation  frequency  and  amplitude,  are  completely  neglected.  This 
l Incarl / at i on  is  probably  valid  for  sufficiently  high  reduced  frequencies,  of 
course,  but  the  extent  to  which  ihts  validity  holds  for  large  deflection 
amplitudes  is  not  clear.  In  the  second  case  the  neglect  of  high  frequency  terms, 
while  self-consistent  for  small  reduced  frequencies,  renders  the  treatment  of 
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last  oscillations  and  gust-like  responses  impossible.  The  confidence  we 
attach  to  low- frequency  results,  not  to  mention  the  exact  consequences  of 
small  disturbance  theory  without  restrictions  to  the  degree  of  unsteadiness, 
are  of  obvious  engineering  interest. 

In  this  paper  two  numerical  methods  are  devised,  the  first  extending 
the  work  of  [2}  to  account  for  the  "back-interaction"  of  the  primary  harmonics 
on  the  nonharmonic  mean  flow,  and  the  second  extending  the  ADI  met  hod  of  [d J 
to  handle  general  unsteady  motions.  These  algorithms,  used  here  to  resolve 
the  issues  formerly  cited,  e.g.,  the  linearity  of  the  unsteady  load  response 
to  changes  in  surface  deflection,  the  influence  of  frequency  on  shock 
excursion  and  mean  shock  location,  etc.,  aim  at  increasing  our  present 
understanding  of  unsteady  nonlinear  interactions.  Apart  from  their 
fundamental  interest,  these  problems  are  important  in  engineering. 

However  we  deal  only  with  oscillating  trailing  edge  flaps  as 
opposed  to  pitching  motions  of  rigid  airfoils.  We  stress  that  these  examples 
are  physically  distinct  and  qualitatively  different.  In  the  former  case  the 
mean  flow  is  roughly  fixed  throughout  the  unsteady  motion  in  the  sense  that 
supersonic  zones  terminated  by  shock  waves  appear  on  both  upper  and  lower 
surfaces.  In  contrast,  for  pitching  oscillations  at  moderate  angles,  each 
ohordwise  surface  alternates  between  shockfree  suberitical  and  supercritical 
shocked  flow  states.  Thus,  to  a certain  degree,  flap  Induced  effects  are 
higher-order.  Qualitative  results  and  trends  obtained  here  for  flapping 
oscillations  generally  do  not  apply  to  pitching  motions,  but  they  may  be 
typical  of  flapping  motions  in  general  and  therefore  of  useful  design  interest. 
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III  ANALYSIS 

Time-Marching  ADI  Method 

Accurate  solutions  for  instantaneous  shock  strength  and  position  are 
needed  to  calculate  net  forces  and  moments,  since  the  departure  of  the  shock 
from  its  mean  position  produces  an  additional  local  lift  contribution 
proportional  to  the  strength  of  the  steady  shockwave  and  the  differential 
upper  and  lower  surface  shock  motion.  These  are  also  essential  to  the 
analysis  of  unsteady  viscid-inviscid  interactions  and  trunsonic  aileron  buzz. 
We  therefore  formulate  the  numerical  problem  for  general  unsteady  motions, 
however,  within  the  framework  of  transonic  small-disturbance  theory.  Consider 
the  nondimen  si  onal  equation  ^A(Pzz  + ~ 4 0^  j)  where  (p  is 

the  disturbance  velocity  potential  normalized  by  c A , c is  the 

ou 

chord,  is  the  freestream  speed  at  infinity  and  is  the  maximum 

thickness  to  chord  ratio.  Let  CO  be  the  oscillation  frequency  and  define  a 

reduced  frequency  by  k = Cd  c/U^  . If  M <*,  is  the  freestream  Mach  number, 

•’2  - 2/3  2 -2/3 

the  required  coefficients  are  A = k“M{1<>  / Sg  , B = kM  M and 

C = (1-M^)/  ^ “S'*  - ( J-  +l)M^c^|  , where  $ = 1.4  is  the  ratio  of  specific 

heats,  nondimensional  variables  having  been  assumed  in  the  form  Z = 11^  kt/c, 
, r 1/3 

i = x/c  and  - da  y/c,  x»  y a*td  1 being  streamwise,  transverse  and 

time  coordinates.  Consistent  boundary  conditions  are  the  specification 

*.  t + 

of  0^  ( l ,0-,  Z ) = f«  4akfj.  on  the  chord  0 - £ 6 1,  where  f is  the 

instantaneous  airfoil  displacement  normalized  by  $ , pressure  continuity 

across  the  trailing  edge  wake  (the  pressure  coefficient  is  evaluated  from 
* 2/3  . . 

C = -2  4 ak ))  and  vanishing  disturbance  velocities  at  inlinity. 

Note  that  a = 0 in  the  low-frequency  approximation  and  a = 1 in  the  general 

case. 
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Solutions  to  the  governing  equation  are  obtained  by  an  alternating- 
dlrectlon-implicit  scheme  which  advances  the  solution  for  <p  at  each 
meshpolnt  from  time  level  Z to  T>  ^ through  the  sequential  procedure 

l - ...op:  ||  J + 

„ . - i wr-f) 


, n?  2 , n 2/3 

where  g = J(C  <p^  -t(l-M— ) ^ / 5Q  )•  In  the  low-frequency  limit  our 

decomposition,  dubbed  GTRAN2,  reduces  to  the  I.THAN2  method  of  £3  J . This 

2 

latter  method  is  0(  ) accurate  and  unconditionally  stable  on  a 

von  Neumann  basis.  For  general  unsteady  motions  the  present  scheme 
consistently  approximates  the  given  equation  with  0(4T)  accuracy;  while 
not  unconditionally  stable,  it  appears  to  be  reliable  and  instabilities 
have  not  been  observed.  Here  $ , <£  and  D are  the  central,  backward 

"V  l i 

and  mixed  difference  operators  defined  in  [3]  ; for  brevity  details  related 
to  tangency,  wake,  shock  und  farfield  conditions,  being  straightforward 
modifications  to  LTRAN2,  will  not  be  discussed. 


N onl i near  Harmonic  Method 

For  small  amplitude  oscillatory  airfoil  motions  the  harmonic  method 
of  W is  adequate.  However,  for  larger  surface  deflections,  nonlinear 
frequency  expansions  must  be  used  to  account  for  the  back-interaction  of  the 
primary  harmonics  on  the*  mean  flow.  On  physical  grounds  the  "mean"  or,  more 
precisely,  the  nonharmonic  part  of  the  totul  disturbance  flow  "freezes"  for 
sufficiently  fast  oscillations  £lj  . Thus  the  dominant  effects  of  nonlinear 
feedback  are  felt  only  for  low-frequency  motions,  where  large  shock 


excursions  are  also  anticipated,  and  we  consistently  set  a-0.  For  unsteady 
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transonic  flows  t he  basic  approach  to  modeling  this  back-interaction  was 
first  described  by  the  present  author  in  [^4^.  Essentially  the  disturbance 
potential  (and  likewise,  the  unsteady  airfoil  displacement)  is  expanded  in 
the  form  P ( | , If  ,Z  ) = p { ,7)  ) +1  t ( ^e1^  ♦ 0^e  iV  ) + . . . where 

higher  harmonic  terms  are  not  shown,  P ^ is  the  complex  conjugate  of  i > *1  ) • 

and  £ is  proportional  to  the  oscillation  amplitude,  and  substituted  into  the 

governing  equations.  For  small  amplitude  motions  a sequence  of  problems  can 

be  obtained  by  equating  coefficients  of  like  powers  of  E ; this  results  in  the 

usual  steady  nonlinear  formulation  for  Pq  and  a linearized  problem  for  j6  ^ 

(e.g.,  see  £‘2  J ).  As  previously  discussed  this  expansion  is  not  uniformly 

valid  in  time:  the  mean  flow  is  artificially  frozen  because  0^  does  not 

depend  on  k or  g . For  slightly  larger  C this  usage  degenerates  even  more 

and  fallacious  results  are  possible.  This  defect  is  remedied  by  equating, 

instead,  coefficients  of  like  harmonics,  in  which  case  an  infinite  system  of 

nonlinearly  coupled  equations  for  the  0^  is  obtained.  To  leading  order 

the  effects  of  secondary  and  tertiary  harmonics  can  be  neglected.  Thus  the 

resulting  equation  for  <fi  ^ is  formally  identical  with  the  equation  obtained 

by  straightforward  linearization;  however,  the  0^  coefficients  are  formally 

unknown  because  the  equation  for  the  "mean  flow"  <t>  now  contains  an 

l T° 

additional  term  of  the  form  aMj/  /d|  which  furnishes  the  required  nonlinear 

2 

coupling.  Although  this  term  is  formally  0(  £ ) it  may  exert  a cumulative 
0(  E ) influence,  especially  near  shockwuves.  The  resulting  equations  for 
P o and  are  solved  iteratively,  in  each  case  using  the  mixed-differencing 
column  relaxation  described  in  [aj  . A trivial  solution  first  initializes 
the  calculations.  One  computational  "sweep"  of  the  flowfield  solves  the  0^ 
equation  subject  to  mean  boundary  conditions  and  an  assumed  zero  0 . Next 


» 
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IV 


the  y ^ equation  is  solved  with  unsteady  boundary  conditions  and  latest 
available  K values,  followed  by  a solution  of  p ^ using  latest 
values,  and  so  on.  Apart  from  the  noted  change  to  the  pQ  equation,  details 
related  to  body  and  wake  conditions,  shock  capture  and  farfield  updates 
follow  £2]  and  will  not  be  rediscussed  here.  We  indicate  the  slow  convergence 
of  the  method,  however. 

RESULTS 

Time-marching  solutions  for  an  unpitched  symmetric  NACA  64A010  airfoil 
section  at  Mach  0.82  were  made  for  an  aft  quarter-chord  flap  with  the 
angular  deflection  history  0<(77)  = l°sin  T at  reduced  frequencies  k = 0.05, 

0.5  and  5.0.  All  time  integrations  were  initialized  with  a steady  flow  and 
marched  over  four  periods  to  rid  solutions  of  undesired  initial  transients. 

gtran2  and  ltran2  results  were  obtained  using  identical  fine  spatial  meshes,  and 
time  steps  were  chosen  in  each  case  such  that  truncation  terms  were  formally 
of  the  same  order.  Computed  unsteady  shockwave  positions  appear  in  Figure  1. 

The  shock  excursions  obtained  in  both  GRTAN2  and  LTRAN2  modes 
for  the  lower  surface  are  similarly  ainusoidal  (second  harmonic  effects  are 
not  visually  perceptible).  Both  suggest  a mean  shock  location  at  56%  chord 
independent  of  frequency  and  recover,  in  agreement  with  Landahl  £l J , the 
shock  motion  freeze  that  obtains  in  the  limit  of  large  reduced  frequency. 
Moreover  GTRAN2  and  LTRAN2  results  are  not  significantly  different,  although 
the  latter  overpredicts  the  peak  excursion  amplitudes  somewhat.  In  the 
k = 0.05  case  the  unsteady  shock  traverses  a net  distance  of  18%  chord,  a 
situation  that  may  be  critical  to  acceptable  boundary  layer  behaviour.  The 
corresponding  net  lift  and  moment  coefficients,  shown  in  Figure  2,  of  course, 
are  similarly  large  in  comparison  to  the  k = 0.5  and  5.0  results.  Figure  2 


■ - 
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also  compares  exact  small  disturbance  results  (GTRAN2)  with  those  obtained 
in  the  low-frequency  approximation;  for  sufficiently  large  reduced  frequencies 


the  latter  yield  inaccurate  amplitudes  as  well  as  phase  information. 

The  chordwise  dependence  of  unsteady  loading  on  flap  deflection  is  best 
represented  using  nonlinear  harmonic  expansions.  Calculations  for  the 
NACA  64A006  airfoil  with  an  aft-mounted  quarter  chord  flap  were  carried  out 
for  Mbo=  0.85  and  kgc  = 0.24  and  Mo.=  0.80  and  kgc  = 0.064  assuming  maximum 
deflection  angles  8 of  1°,  2°  and  3 in  each  case  (kgc  ls  the  reduced 

o 

frequency  based  on  semichord).  We  stress  that  deflections  greater  than  1 , 
approximately,  are  on  the  order  of  the  thickness  ratio  80  . thus  falling 
outside  the  range  where  Ehlers'  straightforward  linearization  ^2^  holds. 

While  our  neglect  of  higher  harmonics  also  implicitly  assumes  "small"  amplitudes, 
the  account  of  mean  flow  distortion  extends  somewhat  the  applicability  of  the 

harmonic  method  to  larger  unsteady  perturbations.  For  convenience  introduce 

r i r i 

the  functions  C , C and  Cn  through  C = C + £(C„  cos  ~C  -C_  sin  T ) and 

Pq  Pi  pl  P P0  P1  pl 


,r.i 


r,i 


define  ^ Cp^  . Figure  3 compares  results  of  Mc 


0.85  and  k =0.24 
sc 


calculations  with  experiment  (^5j  and  with  linear  theory  for  a flap  deflection 

of  £ = 1°.  The  qualitative  agreement  obtained  for  this  supercritical  flow  is 

good  but  precise  results  are  difficult  owing  to  wind  tunnel  wall  interference 

Co  ° 

and  unsteady  viscous  interactions.  Results  of  o = 2 and  3 appear  in 

Figure  4.  The  unsteady  loading  is  seen  to  vary  linearly  with  changes  in 

deflection  angle,  an  observation  not  initially  anticipated  because,  as  previously 

indicated,  § ^ §a  . This  linearity  also  follows  as  a consequence  of  the 

r i 

extremely  weak  dependence  of  C 1 on  £.  Similar  calculations  and  results  for 

P1 

the  Moo=  0.80  and  kgc  = 0.064  case  are  shown  in  Figure  5. 
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V DISCUSSION  AND  CASING  HKMAHKS 

f This  paper  extends  the  scope  of  nil  svsllsble  t lme -marching  Al)l  method  ] 

to  handle  genernl  unsteady  motions  and  also  generallr.es  an  existing  harmonic 
analysis  tool  to  handle  large  amplitude  oscillations.  Immediate  applications 

t using  the  former  method  Include  the  modeling  of  sudden  gust  loadings,  for 

example,  while  t ho  latter  scheme,  for  Instance,  allows  us  to  define  flutter 
boundaries  using  both  frequency  and  amplitude  parameters.  Here  the  new 
® computational  algorithms  are  used  to  study  the  flow  Induced  by  the  transonic 

oscillating  flap  systematically.  The  generalised  ADI  method  reproduces  the 
shock  motion  freer©  In  the  limit  of  large  reduced  frequency,  as  anticipated  by 
^ Lundahl  £l J , and  also  Illustrates  the  slrable  net  lifts  and  moments  that 

ohu rue ter ire  aerodynamical ly  undesirable  low-frequency  osclllutlons.  Moreover 
the  results  of  the  tlme-marchlng  algorithm  developed  here  Illustrate  various 
inaccuracies  associated  with  low- frequency  models  and  suggest  that  their  use 
be  limited  to  reduced  frequencies  (based  on  chord)  of  less  than  0.2,  say. 
Calculated  results  for  the  ^ ■*  1°  case  (e.g.,  see  Figure  1)  Indicate  that 
the  mean  flow  does  not  change  with  frequency. 

The  nonlinear  hurmonic  method  was  used  to  study  the  effect  of  flap 
deflection  on  unsteudy  chordwlse  loading  for  several  low  reduced  frequencies. 

The  inviscid  results,  unhampered  by  the  effects  of  wind  tunnel  wall 
interference  and  unsteady  viscous  Interactions,  Indicate  that  the  local 
unsteady  loading  varies  linearly  with  respect  to  changes  In  deflection 

| 

angle,  even  though  the  deflections  are  not  small.  This  question  was  also 
partially  pursued  by  TIJdeman  . Kxperlments  were  conducted  to  establish 

the  range  of  possible  maximum  flap  deflection  amplitudes  over  which  a linear 
relationship  existed  between  8 and  the  resulting  unsteady  loading. 


l.inoarity  was  est  abl  1 shed  but,  unfort  unately , limitation*  In  tb»>  laboratory 
sot-up  allowed  i>nly  maximum  values  of  1°  for  considerst Ion . For  such  small 
do fleet  Ions  tbe  buffer  furnished  by  the  boundary  layer  reduces  flap 
effectivity  and  the  extent  of  "lnviscld  nonlinearity"  could  not  be  truly 
assessed.  However  the  results  of  the  present  computations  support 
Tijdeman’s  conclusions:  the  observed  linearity  Is  significant  and  may  be  of 
some  consequence  in  the  engineering  of  aircraft  controls.  Tijdeman  further 
reports  that  "the  mean  steady  pressure  distributions  correspond  reasonably 

well  with  the  steady  pressure  distributions  obtained  on  the 

nonoscillating  mixlel"  (this  observation  is  based  on  a flap  deflection  of  l 
with  low  reduced  frequencies  and  supercritical  Mach  numbers).  Again  our 
eomputat ions  support  this  conclusion  (e.g.,  see  Figures  1,  •!  and  It) . 
Additional  numerical  experiments,  not  reported  here,  indicate  that  for  mildly 
supercritical  flows  the  back  interaction  effect  is  a weak  one  and  that 
integrated  normal  forces  and  pitching  moments  respond  sinusoidally. 

The  two  numerical  algorithms  presented  here  for  unsteady  transonic  flow 
with  shockwaves  provide  useful  vehicles  for  unsteady  aerodynamic  calculations 
and  require  only  simple  modifications  to  existing  well  used  computer  codes. 

In  this  paper  they  are  systematically  applied  to  the  study  of  effects  Induced 
by  changes  in  flap  angle  and  oscillation  frequency  in  the  case  of  a symmetric 
unpitched  airfoil.  The  results  agree  with  some  physical  arguments  ot 
l.andahl  jl  ] and  the  experimental  work  of  Tijdeman  [ ft  ] but  they  do  not  bear 
tlie  limiting  restriction  to  small  flap  deflections. 
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FIGURE  3.  MACH  0.85,  REDUCED  FREQUENCY  0.24,  FLAP  DEFLECTION  1 


FIGURE  4.  MACH  0.85,  REDUCED  FREQUENCY  0.24,  FLAP  DEFLECTIONS  1°,  2°  AND  3°  CASES. 
(NONLINEAR  HARMONIC  METHOD) 
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